function f = func(x)
    global G b;
    f = 0.5 * (x') * G * x + (b') * x;
end

